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Abstract 

We investigate ground state configurations of atomic systems in two di- 
mensions interacting via short range pair potentials. As the number of par- 
ticles tends to infinity, we show that low-energy configurations converge to a 
macroscopic cluster of finite surface area and constant density, the latter be- 
• ing given by the density of atoms per unit volume in the triangular lattice. In 

d , the special case of the Heitmann-Radin sticky disc potential and exact ground 

states, we show that the macroscopic cluster has a (unique) Wulff shape. This 
is done by showing that the atomistic energy, after subtracting off a bulk part 
and re-scaling. Gamma-converges to a macroscopic anisotropic surface energy. 

> 

: 1 Introduction 

p 

Q\ ■ The question why a large number of atoms, under many conditions, assembles into a 

^ . periodic, crystalline structure of special geometric shape remains poorly understood 
Q ! from a microscopic point of view. Ultimately this should be a derivable consequence 
of the underlying interatomic interactions, governed by the laws of quantum me- 
chanics. 

In this paper, we address this question in the case of zero temperature, simpli- 
fied interatomic interactions (namely short range pair potentials), and two instead of 
three dimensions. We believe that our methods will remain useful in more general 
settings, including three-dimensional problems. In particular, the prototype case 
allows to understand the key aspects of (1) formation of a local lattice structure (2) 
emergence, as the number N of particles gets large, of a well defined anisotropic "sur- 
face energy" contribution of order 0{N^^^'^~^^) to the atomistic energy of an (exact 
or approximate) atomistic minimizer (3) emergence of an overall geometric (Wulff) 
shape of the atomistic minimizer as a consequence of surface energy minimization. 

We now state our results precisely. Consider particles in two dimensions, with 
positions Xi, ...,xn G M^, whose interaction energy is of pair potential form, 

E{xi, . . . ,xn) = '^V{\xi - Xj\). (1) 



X 
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We will assume that the atomic interaction potential is minimal when |xj — Xj\ = 1 
and short range; see conditions (HI), (H2) and (H3) in the next section for a precise 
statement. One then expects, at least when the potential well of V is sufficiently 
deep and narrow, the following two phenomena to occur: 

(i) E has crystallized ground states for any A^. 

(ii) In the limit N — >■ oo, ground state configurations should assume, in a sense 
to be clarified, a particular overall geometric shape. 

By (i) we mean the following: 

Definition We say that an energy E : (M^)^ M has crystallized ground states 
if its infimum is attained and any minimizer - after translation and rotation - is a 
subset of the triangular lattice 



The energy E given by ([T]) is rigorously known fiUl \TT\ to have crystallized ground 
states when V is given by the Heitmann-Radin 'sticky disc' potential 



or the 'soft disc' potential discussed in [TT]. For analogous results for the closely 
related sphere-packing problem in two dimensions see P El [Hj , and for insights into 
more general potentials see [I3l [7] . 

The present paper focuses on question (ii), which has remained open even in 
the case of energies such as (E]). To describe the overall shape of optimal subsets 
mathematically, we use the following strategy. 

(1) associate to any atomic configuration {xi, . . . ,xn} its empirical measure 



(2) re-scale it so that the total mass and the expected diameter of the support 
of a minimizing configuration remain of order one as — > oo, 

(3) pass to the limit N ^ oc, 

(4) prove that the limit measure is a constant multiple of a characteristic function 
of a set of finite perimeter (i.e., a characteristic function belonging to the space 
BV{R'^)), the constant being given by the density of atoms per unit volume in the 
triangular lattice 

(5) derive, from atomistic energy minimization, a continuum Wulff-Herring type 
variational principle for the shape, via Gamma-convergence, 

(6) use exact solubility of Wulff-Herring variational principles (cf. Taylor [T2] . 
Fonseca-Miiller [8]) to identify the shape. 

The idea not to try and parametrize an atomistic configuration {xi, ..jX^} by 
displacements from a reference configuration (Lagrangian viewpoint), but to study 
its empirical measure YliLi (Eulerian viewpoint) and identify the limit measure 
via a Gamma-convergence result, was recently introduced into the study of many- 
particle energies by Capet and Friesecke in the context of Coulomb systems. 
Note that passage to the empirical measure eliminates the indefiniteness of E under 
particle relabelling. 






+00, < r < 1 
-1, r=l 
0, r > 1, 



(3) 
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The first part of our program, (l)-(4), lias been carried out witliout needing to 
know wlietlier E lias crystallized ground states (in fact, in doing so one needs to 
prove a local form of crystallization). The last two steps, (5)-(6), have been achieved 
assuming crystallization. This leads to the following two theorems. 

The first result is not limited to minimizers, but applies to arbitrary states whose 
energy difference from the ground state is of order 0(A^^/^). 

Theorem 1.1. (Formation of clusters with constant density and finite perimeter) 
Suppose the energy E is given by (QP and the interatomic potential satisfies (HI), 
(H2), (H3). Let {x^^'*, . . . , x^-*} he any sequence of connected (see Definition 3.1 
below) N -particle configurations satisfying an energy bound of form 

E{{xf\...,x^^^})<-6N + CN'/' 

for some constant C independent of N. Let {fJ^N} be the associated sequence of 
re-scaled empirical measures 

N 



i=l 

Then: 

(i) Up to translation (that is to say, up to replacing fijsf by +ajv) for some aat G 
M^j and passage to a subsequence, converges weak* in M.{M?) to jj, E 7W(M^). 

(ii) The limit measure is of the form 

= PXe, 

where p = 2/^/3 (i.e., the density of atoms per unit volume of the triangular lattice 
C) and E is a set of finite perimeter of volume 1/p. 

In fact, any set E of finite perimeter and volume 1/p can occur in the limit, as 
we prove in Section 5. 

Note also that on the atomistic level, quite irregular configurations are admitted 
by our hypotheses. For instance the approximating atomistic configurations may 
contain elastic deformations, cracks and vacancies, or inclusions of phases with dif- 
ferent lattice structure, as long as these only occupy regions of lengthscales smaller 
than A^^/^, see Figure [TJ Theorem II . II say s that on the macro-scale, there neverthe- 
less result well defined clusters of constant, crystalline density. 

For exact minimizers the ensuing cluster has a unique shape: 

Theorem 1.2. Suppose the energy E is given by (J\) and the interatomic potential 
satisfies (HI), (H2), (H3). Assume in addition that E has crystallized ground states 
(as is rigorously known e.g. when V is given by ^). Let , . . . , xj^'' } be any 
minimizing N -particle configuration of E, and let fiN be the associated re-scaled 
empirical measure As N —* oo, up to translation and rotation (that is to say, 
up to replacing p^ by pn{Rn ■ +OAr) for some rotation Rp^ G 5*0(2) and some 
translation vector E M^^) Pn converges weak star to the limit measure 

where h is the regular hexagon conv {±ei, ±62, ±(e2 — Ci)}. 
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Figure 1: Different atomistic configurations satisfying tlie assumptions of Tlim jl.ll 



Figure 2: A ground state of (P), ([3]) for = 72. 

Finally we remark tliat macroscopic uniqueness of tlie limit shape in Theorem II .21 
contrasts with an unexpectedly large amount of non-uniqueness of the discrete min- 
imizers. In a companion paper, we prove that the optimal bound on the difference 
between two ground state configurations and /i'^, 

min{||/i'^ -/i7v(i?- +a)|| : R G 0(2), a G M^} 

(in suitable norms) scales like A^~^/^. Here optimal means that there exists a se- 
quence Nj — > oo of particle numbers for which there is a matching lower bound. Note 
that simple rearrangements of surface atoms only lead to differences of order 2 in 
the associated empirical measures. Hence the A^~^/'^ law shows unexpectedly strong 
fluctuations of finite ground state configurations about the limiting Wulff shape. 



2 Atomistic energy 

Our object of study are low-energy states of many-particle potential energy func- 
tionals of the form 

E{xi,...,xn) = '^V{\xi-Xj\) (6) 

on (M^)^, where Xi, . . . , xat G are the particle positions and the potential V is 
assumed to satisfy the following hypotheses (see Figure [3]): 



V 






-1- 




— ^ — 



Figure 3: Soft Heitmann-Radin potential. 
(HI) (minimum at r = 1) V{1) = —1, V{r) > —1 for all r 7^ 1 

(H2) (behavior at short and long range) There exist constants a G (0, 1], /3 G [1, 00) 
such that V{r) = +00 for r < a, V{r) =0 ior r > P, V continuous on (a,/3) 

(H3) (narrow potential well) The constants a, (3 from (H2) satisfy the condition 
that the ball of radius (3 contains at most six points whose distance from the 
center and mutual distance is > a. 

Simple geometric considerations show that (H3) is always satisfied when a = 1 — e, 
/3 = 1 + e, and e > sufficiently small. The Heitmann-Radin sticky potential ([3]) is 
contained special case [a = (3 = 1). 

The above hypotheses are not aimed at maximum generality, but at simplicity 
of proofs. An important feature is that unlike ([3]), they allow elastic deformations. 

As a simple consequence of the above hypotheses, one has the following lower 
bound on the ground state energy: 

inf E(xi,...,xjv) > -6Ar. (7) 

xi,...,xi^&"S? 

This is because by (H2), (H3), each particle can have a negative interaction energy 
with at most 6 other particles, and by (HI), the interaction energy with each of 
them is > —1. 

Simple trial configurations show that the true ground state energy differs from 
the lower bound by at most 0{N^/'^). More precisely we claim that 

inf E{xi,...,xn) <-QN + aN^''^ + b (8) 

x\,...,xi^&'M? 
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where we may take a = 4-\/3, b = 12. Indeed, consider trial states which are 
hexagonal subsets of the triangular lattice: 



{xi, . . . , xn} = f2 n £, £ as in (I2l), ha C Q C hji, 
where is the closed hexagon of radius R> with "bottom edge" parallel to ei, 

hpi = conv{±i?ei, ±Re2, ±-R(e2 — ei)}, (9) 

is its interior, and R and Q are chosen suitably such that (1 C = N. To infer 
(IH]), it suffices to estimate the energy of these states from above by the right hand 
side of ([8]). For completeness we include the (elementary) argument. 

First, consider the case when Q = h^. In this case the number of "surface 
atoms", i.e. the number of atoms in dhu, equals the length of dhu, i.e. 

Ns = 6R. (10) 

The total number of atoms is obtained by summation over the number of atoms in 
dhr for r < R, 

R 

N = l + ^6r = 3R^ + 3R+1. (11) 

r=l 



Solving equations (ITOl) . (fTTl) for A^^ in terms of yields A^^ = -\/12A^ — 3 — 3. To infer 
the energy, we only need to count the number of neighbors of each atom, where x is 
called a neighbor of y if |a; — y| = 1. Since the 6 "corner" atoms have three missing 
neighbors, and the remaining A^^ — 6 surface atoms have two missing neighbors, the 
energy is 



E{xi, . . . , xtv) = -6A^ + 6 ■ 3 + (A^, - 6) ■ 2 = -6A^ + 2Vl2A^ - 3. (12) 

Now add a partial layer of k atoms in dh^j^i (see Figure ?), where k ranges from 1 
to 6(-R + 1) — 1. Denoting the positions of the layer atoms by xat+i, . . . , x^+k, the 
energy of the configuration is 

E{xi,...,XN,XN+l,---,XN+k) = E{xi,...,Xn) (13) 

k 

+ ^(^-2 ■# (neighbors of XAr+£ in /i/{) (14) 



-1 ■ ^(neighbors of xg in dh^j^i) ). (15) 



1=1 



Here the factor 2 appears because interactions between hn and the layer dh^^i 
appear only once in the sum over i, while the layer-layer interactions appear twice. 
If /c = 1, it is immediate from (|T^ that E{xi, . . . , xat, xat+i) is bounded from above 
by the right hand side of ([H]), so let us assume k > 2. We may arrange the layer so 
that there are at most 5 "corner" atoms with only 3 neighbors, and 2 "end" atoms 
with only 3 neighbors, while all other atoms in the layer have 4 neighbors (see Figure 
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The 5 comer atoms have 1 neighbor in hji and 2 in dh^^i, the end atoms have 

2 neighbors in ha and 1 in dhn+i, and the remaining layer atoms have 2 neighbors 
in each set. Consequently by (fT5l) and ( |T2l) 

E{xi,...,XN,XN+i,---,XN+k) < -E(a;i, . . . ,XAr) + 5(-2 ■ 1 - 1 ■ 2) 

+2(-2 ■ 2 - 1 • 1) + (A; - 7)(-2 ■2-1-2) 
= -6(iV + k)+ I^JVIN -3 + 12. 

Estimating the square root trivially by ^JV2IJ^^P^ yields the desired upper bound 
by the right hand side of ([H]). Although this was not needed here, we remark that 
by the results of [TU] this trial configuration actually forms a ground state of the 
potential (E]). 

3 Compactness and mass conservation 

The results in this section are not limited to minimizers, but apply to arbitrary 
states in which the energy difference from the ground state is of order 0{N^I'^\ See 
Figure [H 

Also, we confine ourselves to connected atomic configurations (see the Definition 
below). In case of disconnected configurations, our analysis can be applied separately 
to the connected components. Note also that minimizers must always be connected. 

Definition 3.1 A finite set 5 C of particle positions is called connected if for 
any two x, y E S there exist xq, . . . ,xn € S such that xq = x, xn = y, and the 
distance between successive points xj-i, xj lies within the interaction range of the 
potential, I < for all j = 1, . . . , A^. 

Also, in the sequel we use the following standard notation. Co(ffi^) denotes the space 
of continuous functions on such that /(x) — > as |x| — > oo, J\A{M?) denotes the 
space of Radon measures on of finite mass (recall M. is the dual of Cq), and a 
sequence of Radon measures /iat is said to converge weak* to /i, notation: fi^ ^ fi, 
if ^2 / dfXN 1^2 f d^i for all / e Cq{^^). 

Proposition 3.1. (Compactness and mass conservation) Assume that the inter- 
atomic potential satisfies (HI), (H2), (H3). Let {x'^^\ . . . he any sequence of 
connected N -particle configurations satisfying the energy hound 

E{{x^l'\...,x'^i;^])<-QN + CN^'^ 

for some constant C independent of N . Let {^n} be the associated sequence of 
Radon measures 1^. Then up to translation (that is to say, up to replacing fij^ 
by /iAr(- + Oat) for some G M.^) there exists a subsequence converging weak* in 
7W(M^) to E A4(M.^). Moreover the limit measure satisfies fiN > 0, J^2 dfi = 1. 

Proof. Throughout the proof, we write S = {x[^\ . . . ,x^^^}, and denote by C a 
constant independent of A^ whose value may change from line to line. 

Since the are nonnegative and have mass 1, they are bounded in A/1(M^) and 
hence, by the Banach-Alaoglu theorem, there exists a weak* convergent subsequence. 



7 



Clearly the limit fi is nonnegative. It remains to show the only really nontrivial 
assertion above, namely mass conservation J^2 dfi = 1. 
The key is to show that 

diamtS = max |a; — ?/| < CN^^'^. (16) 

If such a bound holds, then after translation there exists a fixed ball of radius R 
such that supp/iAT C i?R for all A^. Choose an increasing sequence of functions 
(pn £ (^0(1^^) such that < 0„ < 1, 0„ = 1 on BnR- By dominated convergence of 
(pn to 1, the weak* convergence of /i at to /i, and the fact that supp C Br, 



d/i = lim / (f)ndfi= lim ( lim / (j)ndfij\f) = lim 1 = 1, 

n— too / n— too N^oo J n— too 

completing the proof of the proposition. 

It remains to establish the bound f|T6|) . We begin by introducing a notion of 
neighbors and a notion of local energy. Let e G (0, 1), small enough so that condition 
(H3) is satisfied for a = 1 — e and (3 = 1 + e. 

Definition We say that x G Sat is a neighbor of y G S'tv if |x — y| G [1 — e, 1 + e]. 
The set of neighbors of a; G 5 is denoted by M{x). 

Also, independently of the above notion of neighbors, we define a local energy, as 
follows: 

y£S\{x} 

SO that 

E{S) =J2Eioc{x). 

xes 

By hypotheses (H1)-(H3) on the potential V and the finiteness of the energy of Sn, 
it follows that 

#Ar(x)<6, Eeocix)>-Q, (17) 

with both bounds being sharp as is seen by considering points x in the lattice C. A 
key point now is that when the number of neighbors of x is not equal to 6, the local 
energy of x is bounded away from its optimum by a finite amount: 

#M{x) < 6 ^ Eioc{x) > -6 + A, with A := min V{r)-mmV{r) > 0. (18) 

Next, we construct an appropriate set in the plane associated with the configuration 
S. For each x G 5, let V(x) be the Voronoi cell of x, 

V(x) = ^^y em^ : \y - x\ <\y - x'\ for all x' G 5\{a;}|. (19) 

As V{x) may be unbounded, it is useful to introduce in addition its truncation 

Vtruncix) ■.= V{x)nB,{x), (20) 

where Br{x) is the ball {y & M."^ : \y — x\ < 1} . We then define 

^ ■■=[j Vtruncix). (21) 
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Elementary geometric considerations show that when e > is chosen sufficiently 
small, then 

=Q^dnn dVtruncix) = 0, (22) 

that is to say the cells associated to points with the maximum number of neighbors 
do not contribute to the boundary of fl. Finally, since by construction Vtrunc{x) C 
Bi{x) and Vtrunc{x) is convex, we have the following bound on the length of its 
boundary 

\dVtruncix)\<2n. (23) 

Consequently, denoting 

dS:={xeS: #Ar(x) < 6}, (24) 

using the plausible fact proved in Lemma 13.21 below that due to the connectedness 
of S the set Q is connected, and 



diam5 < diamfi < l\dn\ < ^ ^ \dVtrunc\ < vr #95. (25) 

On the other hand, recalling N = and using the assumption on E{S) in the 
lemma and ffTSl). 



- 6iV + ^ ^^^^ ^ ^ ^ ^Q^^ ^26) 

and consequently 

#55 < ^N^'^. (27) 

Combining flTTj) establishes (ITBI) . The proof of the proposition is complete, 

except for the following elementary geometric lemma. □ 

Lemma 3.2. // a configuration S is connected in the sense of the definition in 
SectionlM and the energy E{S) is finite, then the set Q defined in ^21\} is connected. 

Proof. We begin by establishing an elementary inequality relating the constants a 
and P appearing in the hypotheses on the potential V. By (H3), the ball of radius 
/? around the origin does not contain the regular heptagon around the origin with 
sidelength a, or equivalently 

— > p sm — , = — . 
2 ^ 2' ^ 7 

Hence 

/3< ^ . , , = 1.152382... a. (28) 
2 sm0/2 

To prove the lemma, it suffices to show that if |x — ?/| < f3, the truncated Voronoi 
cells around x and y have nonzero intersection. To establish this, it is enough to 
show that the midpoint m = belongs to Bi{x)nBi{y) and V{x)nV{y). The ffist 
inclusion is immediate from |m — x| < /3/2, ( 1281) and a < 1. The second inclusion is 
equivalent to 

I I ^ I II h 

\m — z\ > (= \m — x\ = \m — y\) 
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for all z G S\{x^y\. Now any such z belongs to M^\(5a(3;) U Ba{y) (otherwise the 
energy of the configuration would be infinite), and any closest point p to m in the 
latter set has distance \p — m\ = a'^ — {\x — y|/2)^, so it suffices to show that 
i/a^ — (|a; — y\/2y > \x — y\/2, or equivalently 

« > (29) 

But thanks to ( !28l) we have 

\x-y\< f3 < 1.152382... a < \p2a. 
This establishes ( !29l) . completing the proof of the lemma. □ 



4 Proof of formation of clusters with constant 
density and finite perimeter 

We now prove Theorem II. H i.e. we show that under the assumptions of Proposition 
I3.lt the limit measure /x of the re-scaled empirical measures of atomistic configura- 
tions is a constant multiple of a characteristic function., the constant being given 
by the density of atoms per unit volume of the triangular lattice. 

Proof of Theorem 11.11 As in the proof of Proposition 13.11 we use the decom- 
position of into truncated Voronoi cells Vtmncix) (see (IT9l)- (l20l) ) associated with 
the points x of an atomistic configuration S. The main technical idea of the proof 
is to introduce and investigate the following volume excess function 

^^"^^ '-^\\Vtrunc{x)\~ \h[^^^ix)\[ ^^^^ 

Here h[.{x) is the closed hexagon of radius r around x which is obtained from hr 
(see dH])) by a 30° rotation around the centre. This hexagon is the Voronoi cell of 
any interior point x of a particle configuration S on the triangular lattice i^. 

Besides the sequence of re-scaled empirical measures (jl]) , we will make use of the 
following auxiliary sequences: 



JV lAf-l/2T/; (^M' ^ ' 



N ^ \N-^''Vtrunc{x)\ 
xGSn 



Step 1 First we claim that 



N ^ \N-^/^h' Jx)[ ^ ' 



/ijv -/ijv ^ in A1(M2). ^33^ 
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Indeed, since for any xq G Sjy, Vtrunc{xo) C Bi{xo), we have for all test functions 
e Co(M2) 

1 



N-y^Vtrunc{Xo)\ 

N-^/2Vtrunc(xo) 

x) - (f){y)\ 



(l){xo/vN) - dx 



< sup 

\x-y\<N''^/^ 



and consequently 



[ (l)dfiN- ! (pdfiN < T7 sup \(f){x) - (t>{y)\ 

sup |</.(x)-0(y)| ^0 (iV^oo), 

|a;-j/|<Ar-i/2 

establishing (133!) . 

Step 2 Next, we compare the measures fiN and /lAr. We will show that 

/i^v - /iiv ^ in L^(R^) 



(34) 



(and hence, a forteriori, Jijq — /ijv ^ in ^A{M?)). First, we decompose both 
measures into their "interior" and "boundary" parts, as follows. Here OSn is as 
defined in ( l24l) . and intiSAr := SnXOSn- 



Hn — Aat + un, X 



N 



- E 



and analogously for /i^v. Roughly speaking, we will argue that the difference between 
uisf and i>Ar is small because they are small separately, due to the fact that the number 

of boundary atoms grows only like N^^^, and that the difference between A^v and 
Aat is small because otherwise this would cost elastic energy. 

To make the first argument precise, we use that for any configuration S = 
{xi, . . . ,xn} with finite energy, Vtrunc{xi) 3 Ba/2{xi), due to the fact that |xj — Xj| > 
a for all j i, so that the nearest point in 5 to y G Ba/2{xi) is Xi. Hence 



XN-^nVtrunc{x) 



\B. 



a/2\ 



\h' 



=:K 



Since ^OSn < {C/A)N^/^ (see (EZD) and Vtmndx) C Bi{x), we infer that 



\'^n-'^n\l^ 



< ^iVV2 max ||XAr-i/.y,_(x)IU^ < ^^''^'l^il - (iV 00) 

3; to A' 



CK 



(35) 
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Next we analyze the difference between the interior parts, 

1 1 



Aat — Aat — 



For any point x G intiSAr, #A/'(x) = 6, so 7V(x) = {t/i, . . . , ^/e}, where we may 
assume that the Uj are numbered so that yj — x = rj(cos0j, sin0j), < 0i < 02 < 
... < 06 < 271. For elementary geometric reasons, namely that the only way to 
arrange the i/i so that \yj — x| = 1 for all j and \yj — yj-i\ = 1 for all j is to place 
them at the corners of a regular hexagon around x, we have 

6 

\ytrunc{x)\ ^ \h^/^{x)\ if ^ V {\x - y\) 

yeN{x) j=l 

Hence by continuity, given 6 > there exists A (5) > such that for all x G intiSAr 
the volume excess ( 130|) satisfies the following implication: 

> 6 

6 

^Eeoc{x):= J2 ^(k-l/|) + 5^V^(|%-+i-%-|)>-12 + A. (36) 

yeAf{x) j=l 

It will be convenient to extend the function Eioc{x) to all of Sn, by setting it equal 
to —12 when x G SiSat. Combining the energy bound assumed in Theorem 11.11 the 
fact that each particle pair appears in Egodx) for at most four x ^ Sn whereas it 
appears twice in E, and ( l36l) yields 



2 

xes 



> i (^-12#|x e : X e OSn or ^(x) < 
+ (-12 + A)#|x e int Sn : $(x) > 5} 

> -6A^ + G intcSTv : $(x) > 5|. 



Consequently 



and so 



#|x e int5jv : $(x) > < '^N^/^ 



IAtV - A7v||li(k2) - / I 5Z ^N-^/^Vtruncix)(TTT- 

J I \ I V4rr.^,. 

a;Gint (Sat 



^trnnc(a;)| |/ii/v^(a;)| 



< Yl ^i^)\\XN-y^Vtruu.{x)\\L^ 
xGint cSjv 



xSint5jv 

< N''^\Bi\(6N+ max <l>(x) #|x G int 5 : ^(x) > ^j) . 

V xSintiSjv .1 J / 
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Letting N oo gives 

limSUp IIAat - AAr||Li(R2) < \Bi\S. 

Since 6 > was arbitrary, the left hand side equals zero, that is to say \n — 
in L^(M2). Together with (EHl), this estabhshes (IM]). 

Step 3 Having established that the limits of /iat, /lAr, fi^ all coincide, it suffices to 
study the limiting behavior of the sequence fl^, which - as we shall see - is compact 
in a much "stronger" space. After a change on a set of measure zero, 

P'N = 777 rXjV-i/anjv' = [J Vtruncix). 

' l/ys' x€Sn 

Clearly the fiN only take values in {0, p}, belong to the space are bounded 

in L^(M^), and - by ( |T6l) - after translation are supported in some fixed ball of 
radius R. We claim that they are also bounded in BV{R'^). This is because \dflN\ is 
bounded by a constant times N^^'^, by (12^ and (E7|) . Consequently by the Banach- 
Alaoglu theorem and the compact embedding BV{B^ ^ L^{Bji), a subsequence 
converges weak* in BV and strongly in to some limit fl G BV. By the strong 
convergence, ft only takes values in {0,p}, and by Steps 1 and 2, /i = /i. The proof 
of Theorem 11.11 is complete. □ 



5 Gamma- convergence of surface energy and emer- 
gence of Wulff shape 

In case of low-energy states, we now study the shape of the limiting cluster E 
obtained in the previous section. 

Our current methods are restricted to exact minimizers and to energies with 
crystallized ground states. Here the problem simplifies because by (HI), (H2) and 
(H3) we may without loss of generality assume that the interaction potential is given 
by the Heitmann-Radin potential (E]) (see Figure HI). 



V 



-1 



+00 
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r 









Figure 4: The Heitmann-Radin 'sticky disc' potential. 
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For such configurations, we can derive a hmiting variational principle, as follows. 
It is useful (as in [1]) to re- formulate the minimization problem for the atomistic 
energy E in terms of empirical measures instead of particle configurations. Define 
the following energy functionals on the set V of probability measures on (i.e., 
nonnegative Radon measures on M? of mass 1): 

[ k^v^a, NV{Ny^\x -y\)d^^di^, /i = ^ Ef=i S^./VN 
/Ar(/i) := < for some distinct Xi E C (37) 

[+00, otherwise. 

This definition says that I^ifJ^N) = E{xi, . . . , xn) when /iat is the re-scaled empirical 
measure (jll) of the configuration {xi, . . . ,xn}- In particular, the re-scaled empirical 
measure minimizes In if and only if the underlying configuration minimizes E. 
We now show: 

Theorem 5.1. The sequence of functionals N~^^'^(In + 6N) Gamma-converges, 
with respect to weak* convergence of probability measures, to the limit functional 
loo ■ — > M U {oo} given by 

{jQ,^T{i'E)d'h}{x), /i = -^Xe for some set E of 
finite perimeter and mass (38) 
+00, otherwise, 

where T is the function 

rM^2(.,-±|)fa.^(-j;).,.|o.fl, (39) 

extended -periodically . 




Figure 5: Plot of the limiting surface energy density F in dependence of tp and polar 
plot of F clS cl function of the normal u. 



The physical idea which is made quantitative here that crystalline anisotropy 
gives rise to nontrivial surface energies goes back to Wulff (1901) and Herring (1951) 
in the physics literature. Rigorous interfacial energy results have a large and sophis- 
ticated statistical mechanics literature, particularly in the context of the Ising model 
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(see e.g. the book by Dobrushin, Kotecky, and Shlosman [0]). For a previous treat- 
ment of an interfacial energy problem via Gamma-convergence see [1]. 
Proof. Lower bound. Let Sn = ■ ■ ■ , x^n^} C £ be any sequence of A^-particle 

configurations whose associated sequence of Radon measures fiN (see ([4])) weak*- 
converges to a probabihty measure fi We need to show that hminf JAr(/iAr) > 

/oo(/U)- 

Associate to ^Stv the following auxiliary set 

where h'^{x) denotes the open hexagon around x of radius r introduced in the pre- 
vious section. Note that for x E C, the closure of h'^^^i^) is the Voronoi cell of x 
with respect to the complete lattice C 

The boundary dH^^ is a disjoint union of simple closed polygons Vi, . . . , Vm- 
Because the boundary dH^ oscillates on the atomic scale, we define yet another 
auxiliary set H'^ (see Figure E]) which removes these oscillations and will hence allow 
us to obtain a sharp lower bound on /Ar(/XAr) via standard weak lower semicontinuity 
results on surface functionals. If Vj C dHf^f is a simple closed polygon, then Vj = 
UiliN+i)^j] with segments of length l/^/SN for 

vi, ...,Vm.e ( -^(ei + ea) + -^c]u( -^{2e2 - ei) + -^c] , v^+i := Vi. 



3VN VN J \3VN VN 

As the corner points Vi alternate between the lattices j^i^i + ^2) + '^'^ 
j^{2e2 — ei) + -^jC comprising the dual lattice of £, m is an even number, and 
Vj = [J'^i[v2i-i, V2i+i] is a simple closed polygon. We now define the set H'j^ C 
as the unique closed set with G H'j^ such that dH'j^ = [Jjli It can be easily 
seen that 

< (40) 

Here, for two sets A and B, AAB denotes the symmetric difference {A\B) U {B\A). 
Note that [v, w] fl [x, y] 7^ defines a one-to-one correspondence between boundary 
segments [v, w] of and nearest neighbor bonds [x, y] of -^'C., i-e. x,y E 
with \x — y\ = such that x G -^Sn and y ^ -^Sjy. But then there is a one-to- 
two correspondence between boundary segments [v, w] of H'^^^ and nearest neighbour 
bonds [x, y] such that x E S and y ^ S. Since boundary segments of H'p^ are of 
length -t=, this shows that 

/7v(/iiv) + 6iV = 2^/Nn\dH'j^). (41) 
Precisely as in Step 1 in the proof of Theorem 11.11 we see that 

1 Xjv-i/2fe;^^(a;) 

^= iv ^FW— (42) 
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Figure 6: Constructing H'^^ from H^. 



weak* converges to /i. Now — similarly as in Step 3 in the proof of Theorem 11.11 — we 
obtain that after a change on a set of measure zero 



Moreover, if -^(/^^(/iyv) + 6A^) is bounded it follows from (1411) that the restrictions 
/^Aflsfl are bounded in BV{Bb) for each i? > 0, and hence converge to yu|_B^ in 
L^{Bji), due to the compact embedding BV{Br) > L^{Br). Since ||/i||ii(R2) = 
||/iAr||Li(K2) = 1 for all iV, we obtain that even yU^v yU in L-'^(]R^) and jj, = pxe for a 
set E of finite perimeter and area K 

Furthermore, if (a subsequence of) -^{iNiP'N) + 6A^) is bounded, then by (1^ 
and I^Tij PXh'^ — ^ = PXs in i^H^^)' whereas by fHT]) . flHUl) and the fact that each 
boundary segment of H'j^ is parallel to ei , 62 or ei — 62 

{InM + 6iV) = 2n\dH'^) = j T{v) dn\ 



N 

Now since F, extended to a 1-homogeneous function on M", is convex and satisfies 
a growth condition of the type > c|z/| for some c > 0, a lower semicontinuity 
result for SBV functions (cf. Theorem 5.22 in Ambrosio-Fusco-Pallara [2]) establishes 
the lower bound 



liminf ^(Jjv(/ijv) +6iV) > \ T{uE)dV}\ 

N^oo i/jV J 

d*E 



X] 



Upper hound. The idea to find a recovery sequence is to approximate gradually a 
set of finite perimeter E by sets of simpler geometric shape. To be more precise, 
a set of finite perimeter will be approximated by C°° sets, C°° sets by sets with 
polygonal boundaries and polygonal sets by polygonal sets having all their corners 
in -C for some n G M. The recovery sequence for E will then be extracted through a 
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diagonalization process. Finally, a suitable continuity property of surface integrals, 



lim / T(u(x))dnUx) 

N^oo J ^ ^ ^ ' 

d-PN 



V{u{x))dn\x) 



d*P 



where (-Pat) is the approximation sequence and P denotes the approximated set in 
the respective step, will complete the proof. 

Step 1 First, let P C be a bounded set with polygonal boundary such that every 
corner of dP lies in for some n G N. Assume that the volume of P is ^ + 
Consider the sequence of particle configurations Sn,N = C H \/NP consisting of 
M = Mn^N atoms. Elementary geometric considerations show that 



(43) 



\M^,N-N\<c{anN + ^N). 



for a constant c independent of n and A^. Let ^n.N denote the associated rescaled 
empirical measures. Clearly, /in,Af PXp as iV — > oo. In addition, a straightforward 
calculation of the associated energy by evaluating the surface energy contribution 



N ^ 



X 



along line segments of dP gives 



(/A/(/in,Jv) + 6M) - / V{u{x))dn\x 



dP 



< c a„ + 



(44) 



for a constant c independent of n and A^. Indeed, an elementary argument shows 
that if 5* is a boundary segment of P of length L with normal u, then the number 
of pairs {x,y) G Sn,N x (-^ \ <Sn,N) with \x — y\ = 1 such that the segment [x,y] 
intersects \/N S is equal to 

L 



[VN\T{iy)L + 0{L) = VN { T{u)L + O 



Step 2 Now let P C be any set of volume ^ with polygonal boundary. By 
perturbing the corners of dP slightly, it is easy to see that there is a sequence of 
polygonal sets P„ whose boundary vertices lie in -C satisfying 



|P„AP| — s> as n — s> oo 



(45) 



and 



r(i^p„) J r(z/p) as n ^ oo. (46) 

a'Pn d*P 

Now choosing n = n{N) oo appropriately, we obtain a sequence of configurations 



S 



N 



5, 



n{N),N 
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from the configurations constructed in Step 1 such that the associated rescaled 
empirical measures jiN satisfy fiN ^ PXp aiid 

^ {Im{~Pn)+QM)^ I T{u{x))dn\x) 



N 

dP 

by (USD, dSD and 

Now if ^iSat 7^ A^, then the configuration Sn needs to be modified either by 
adding or removing some set of points with neghgible surface term. By (H3l) -^\M — 
N\ Q as N ^ oo. If Sn contains less than particles, then we can add the missing 
number of points on lattice sites within some parallelogram whose sidelengths are of 
order \/N — M and such that the points in the parallelogram are a distance greater 
than 1 away from S^. Then none of the new points interacts with S^. In case 
M > N, just remove M — N points of Sn lying in a common parallelogram whose 
sidelengths are of order \JM — N . (This is always possible if is large enough.) 

Denoting the corresponding rescaled empirical measure hy and recalling that 
-i:|M — A^l ^ it is not hard to see that fijy ^ pxp N ^ oo. Since furthermore 



ihiiM + 6M) - -y={lNM + 6iV) 



for some constant c > 0, we also obtain 

1 



<c\lj^{M-N) 



dP 

by (USD and (HH). 

Step 3 A recovery sequence for a general set E of finite perimeter is now obtained 
by a diagonalization argument due to the following density and continuity results. 

Suppose first E is a. bounded set of finite perimeter with C°°-boundary and 
volume By piecewise linear approximations of dE and scaling with factors close 
to 1 we easily construct approximations P„ which have the same volume, a polygonal 
boundary and satisfy 



Xp„ ^ Xe and J T{upJ J T{ue). 

a*p„ d*E 

Now let £■ be a bounded set of finite perimeter with volume As shown e.g. in 
[3], Proposition 4.7 and Remark 4.8], there are bounded sets En of finite perimeter 
with C°°-boundary and — after rescaling with a factor close to 1 — volume ^ such 
that 

Xe„ ^ Xe and J T{ueJ J T{ue). 

d*En d'E 

Finally note that a truncation argument yields that an analogous result holds for 
approximating sets of finite perimeter with bounded sets of finite perimeter. This 
concludes the proof. □ 



To complete the proof of Therem II. 2[ it remains to infer convergence of mini- 
mizers. A technical detail we need to pay attention to is that, unlike in many other 
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Gamma-convergence results, here sequences with bounded energy are not in general 
compact in the topology in which the Gamma-convergence occurs. This is because, 
due to the translation invariance of the functional N-^/^{In + 6iV), sequences can 
lose mass as ^ oo. 

Proof of Theorem 11.21 Let {x"-^\ . . . be any minimizing A^-particle con- 

figuration of and let /^at be the associated re-scaled empirical measure (jlj). By 
the connectedness of minimizing configurations and Proposition 13.11 after suitable 
translations /iTv ^ fJ^N{- + ciN) the limit measure /i has full mass. Hence sequences of 
exact minimizers are compact in the topology in which the Gamma-convergence oc- 
curs (namely weak* convergence of probability measures). By standard arguments 
in Gamma- convergence, is a minimizer of the limit functional loo- 

We now appeal to the uniqueness theorem for Herring type energies due to Taylor 
(in the language of geometric measure theory), in the version by Fonseca and Miiller 
(who work in the present setting of boundary integrals for sets of finite perimeter): 

Theorem 5.2. fTIi [S] A functional of form 

I{E) = j T{u{x))dTT-\x), 

with r : S^''^^ [0, oo) continuous and bounded away from zero, is minimized 
over sets E C of finite perimeter and volume 1 if and only if E agrees, up to 
translation and up to a set of measure zero, with XWr, where Wr is the Wulff set 

Wr := \^x eW : x-iy < T{u) for all v G 

and A > zs the unique normalization constant such that AH^r ^os volume 1. 

In the present case of the energy fl39|) . an elementary calculation shows that the 
Wulff set is given by the intersection of the six half-spaces x ■ v < r(z/) for the 
minimizing normals V2-Kj/&i j = 1? • • • ,6, i.e. a regular hexagon with bottom face 
parallel to the ei-axis (see Figure [T]). This completes the proof of Theorem 11.21 




Figure 7: Polar plot oi u ^ r(z/) and the associated Wulff shape Wy- 
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